In-Medium Similarity Renormalization Group for Open-Shell Nuclei 
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We present a nonperturbative derivation of effective valence-shell Hamiltonians in the framework 
of the recently developed in- medium similarity renormalization group (IM-SRG). As a first appli- 
cation, we calculate the spectra of p- and sd-shell nuclei, ^Li and ^^O, based on evolved chiral 
nucleon-nucleon interactions. For ^Li, the spectrum is in very good agreement with ab-initio re- 
sults. For ^^O, the IM-SRG provides a new method for the shell model to systematically go beyond 
effective interaction techniques based on diagrammatic expansions. 
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Introduction.- Advances in ab-initio methods for nu- 
clear structure combined with nuclear forces based on chi- 
ral effective field theory (EFT) have lead to many excit- 
ing developments for light nuclei and medium-mass nuclei 
around closed-shell configurations (see, e.g., Refs. [1]). 
For open-shell systems with many valence nucleons, how- 
ever, the shell model remains the most successful ap- 
proach to understand and predict nuclear structure, in- 
cluding the evolution of shell structure with changing 
neutron and proton numbers, properties of ground and 
excited states, and electroweak transitions [2|. More- 
over, the shell model has recently revealed new insights 
to the impact of long-range tensor [3] and three-nucleon 
forces [J] in neutron-rich nuclei. These are dominated by 
pion exchanges, which provides a link between nuclear 
structure and developments in chiral EFT interactions. 
Despite the many successes of the shell model, the micro- 
scopic derivation of effective interactions and operators 
among valence nucleons from nuclear forces is still largely 
based on perturbative approaches where the convergence 
remains an open problem. 

In this Letter, we present a new nonperturbative 
derivation of effective valence-shell Hamiltonians in the 
framework of the in-medium similarity renormalization 
group (IM-SRG), which we recently developed for closed- 
shell nuclei [5, 6]. The IM-SRG is based on a renormal- 
ization group evolution that decouples degrees of freedom 
that are not relevant for the problem of interest. We show 
how the IM-SRG can be generalized to open-shell systems 
away from doubly- magic nuclei. For ^Li, we present first 
results for the ground-state energy and spectrum in very 
good agreement with ab-initio methods. We then dis- 
cuss for ^^O how the IM-SRG goes significantly beyond 
effective interaction techniques based on diagrammatic 
expansions [7], opening up a promising new method to 
connect nuclear forces and the shell model. 

IM-SRG and generator choices.- The IM-SRG starts 
from a Hamiltonian H that is normal ordered with re- 



spect to a finite-density reference state |^) (e.g., the 
Hartree-Fock ground state): 

H = Eq^Y^ fij {alaj} + ^ XI ^^Jk^ {ala]aiak} , (1) 



ijkl 



where the normal-ordered strings of creation and annihi- 
lation operators obey (<^|{a| ■ ■ ■ aj}\^) = 0. We include 
normal-ordered 0-, 1-, and 2-body operators, Eq, /, and 
r, which approximately include induced 3- and higher- 
body interactions, and solve the IM-SRG flow equations 
to obtain the evolved Hamiltonian H{s) [5]. We refer to 
this truncation as IM-SRG (2) since we keep up to normal- 
ordered 2-body operators. The evolution is equivalent to 
a series of unitary transformations that are designed to 
evolve H{s) diS s ^ oo to an appropriately defined "di- 
agonal" part H'^is) [i,IH: 

H{s) = U{s)HU'^{s) = H'^is) + H'^'^is) -^ H'^ioo) . (2) 

The unitary transformation U{s) is determined by the 
generator 77(5) = [dU{s)/ds] U^{s), which is constructed 
from the diagonal part. 



r?(s) = [iJ^(s),iJ(s)] = [iJ^(s),ff°^(s)] 



(3) 



and guarantees that the "off-diagonal" coupling H'^^ is 
driven to zero with increasing s. 

For the ground state of closed-shell nuclei, one elim- 
inates all terms that couple the reference state |<l>c) to 
the rest of the Hilbert space. This is achieved when the 
matrix elements between |$c) and all n-particle-n-hole 
states vanish, {npnh\H {00) \^c) = 0. Therefore, one 
takes H^^ to be composed of all 1- and 2-body opera- 
tors that connect hole (h) with particle (p) states so that 
{H^^} = {fph^ rpp//i/i/| plus hermitian conjugates, as was 
demonstrated in Ref. [5]. 

For open-shell nuclei, particle states p can either be va- 
lence particles or particle states above the valence space. 
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FIG. 1: Schematic illustration for the valence-space decou- 
pling by the IM-SRG evolution from H{s = 0) ^ H(oo) for 
two valence nucleons. 



which we denote by v and q respectively. We want to de- 
couple states that are not in the valence space, spanned 
by \^y) = {al^ • • -al^} \^c)i where TV is the number of 
valence nucleons and \^c) is the reference state for the 
core nucleus with Ac nucleons. This can be realized by 
defining H"^^ through the set of matrix elements 

l-^l J ^ XJph^ Jqv-) J- pp'hh' -) J- pp' {vh or hv) -) ^ {pqov qp)vv' j -> 

(4) 
where p = v^q^ plus hermitian conjugates. As an alter- 
native generator choice i^2^, we also drive the one-body 
part to diagonal, so that we define 



{^2 } — {^1 ^fpp'^fhh'}' 



(5) 



These generators both lead to a diagonal part Hf 2 where 
states outside the valence space are decoupled by the IM- 
SRG flow, illustrated in Fig. [TJ leading to 



PH^2{^)Q = QH^ 2(00) P = , 



(6) 



with P = ^^\^y) {^y\ said Q = 1 - P. The ofT-diagonal 
parts in Eqs. (j4]) and (J5j) can also be derived using the 
counting operator C = ^iCi{alai}^ with Ci = 1,0,-1 
for q^v^h states, respectively. The C operator counts 
the number of excitations on top of a valence-space 
state \^v)- It is then straightforward to verify that 
the above choices of Hi^ ensure = CH {00) \^y) = 
[C, H{oo)] \^v)^ which leads to the decoupling of valence- 
space states from arbitrary excitations. 

After the IM-SRG (2) evolution, the effective valence- 
shell Hamiltonian is given by Hq^^ = PHf2{^P ~ ^0 ""' 
where Eq "" is the 0-body piece of the evolved Hamiltonian 
corresponding to the ground-state energy of the core. We 
then solve a reduced eigenvalue problem in the N valence- 
particle space. 



iJeff \Xn) = {E^ - Ei 



\Xn 



(7) 



Results.- We next present first applications of the IM- 
SRG to two open-shell nuclei, ^Li and ^^O, consisting 
of two valence nucleons on top of the closed-shell nuclei 
^He and ^^O. All results are based on the SRG-evolved 
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FIG. 2: The ground-state energy of Li versus harmonic- 
oscillator parameter hio obtained by a diagonalization of the 
IM-SRG (2) ifeff in the p-shell, using generators H'f and i7§^ 
(left and center panels). The flow equations are solved in 
the HF basis truncated to Cmax = max(2n + /) single-particle 
excitations. For comparison we show in the right panel the 
convergence with increasing A^max of the NCSM energy jl2 ]. 
where the dot-dashed band is the extrapolated result. 



(A = 2.0 fm-^) N^LO NN potential of Ref. [To^. We 
begin with ^Li, which is a sufficiently light nucleus to al- 
low a direct comparison of our IM-SRG (2) results with 
exact diagonalizations of the 6-body problem using the 
No-Core Shell Model (NCSM). For the IM-SRG calcu- 
lations, the flow equations are solved in the Hartree- 
Fock (HF) basis of the initial Hamiltonian truncated to 
Gmax = max(2n + /) single-particle excitations. The re- 
sulting JT^eff is then diagonalized in the p-shell to obtain 
the ^Li energy levels with respect to the ground-state en- 
ergy of the ^He core. There is a subtlety that arises due 
to the self-bound nature of atomic nuclei. As we wish to 
minimize spurious center-of-mass motion, we work with 
the intrinsic Hamiltonian li\rv\. = H — 'P'^/(2mA) where 
P = ^^ Pi and A = 6 for ^Li. Therefore, the eigenval- 
ues of Hq^^ correspond to the excitation energies of ^Li 
with respect to the ground state of the unphysical "^He 
nucleus obtained using H\^t with A = 6. Consequently, 
to get the absolute ground-state energy of ^Li, we do a 
separate IM-SRG (2) calculation of the ground state of 
the unphysical ^He core and add this to the eigenvalues 
of i^eff • We have checked that for large Cmax spaces, the 
center-of-mass factorizes as in Ref. jllj . 

Figure [2] shows the convergence of the ground-state 
energy of ^Li with increasing emax excitations. The left 
panels give the IM-SRG (2) results using the two differ- 
ent generators H^"^ and i^2^, while the right panel shows 
NCSM energies p^ for comparison. Since the single- 
particle Cmax truncation is different than the NCSM 
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FIG. 3: Convergence as a function of Cmax of the excitation 
energies of ^Li obtained by diagonalizing the IM-SRG(2) ifeff 
in the p-sheh using the Hi^ generator. The HF basis at a 
fixed hu = 24MeV is used for the IM-SRG(2) calculations. 
For comparison we show the convergence with A^max of the 
NCSM energies at the same huu value J12l] . 



^max, the convergence pattern of the two methods is 
expected to be different (an emax space for the same 
value is substantially larger). The IM-SRG(2) ground- 
state energy converges to — 32.7(3)MeV, where contri- 
butions from normal-ordered three-body interactions are 
expected to be repulsive (similar to triples correction in 
coupled-cluster calculations) [5] , in very good agreement 
with the extrapolated NCSM value -32.0(2) MeV [12]. 
The generator dependence in Fig. [2] is found to be very 
weak, indicating that the error from truncating the IM- 
SRG equations to two-body operators is indeed small. 

Next, we study the convergence properties of the low- 
lying excited states of ^Li. The left panel of Fig. [3] shows 
the convergence of the IM-SRG(2) spectrum as a func- 
tion of emax at a fixed value of hu = 24 MeV. How- 
ever, the hu dependence is very weak in the HF basis 
for large emax- The right panel shows the convergence of 
the NCSM spectrum with A^max- The low-lying 3+, 0+ 
and two 2+ states converge rather well and are in reason- 
able agreement with the NCSM results. The high-lying 
1+ state is not yet converged even at the largest space 
for both methods. As the convergence is very poor in 
the NCSM with a harmonic-oscillator basis, this could 
indicate that this state has an extended structure. 

Our results for ^Li are very encouraging and show that 
the IM-SRG provides a new method to derive effective 
valence-shell Hamiltonians that accurately reproduce the 
low- lying spectrum obtained with ab initio methods, but 
at a polynomial scaling ~ N^N'^ with the number of hole 
and particle orbits. Recently, other methods have been 
explored for open-shell nuclei, including two-particle at- 
tached coupled-cluster theory [13,], which however leads 



to nonhermitian effective Hamiltonians, and the NCSM 
with a core [14], which requires a NCSM solution of the 
full problem and is therefore limited to lighter nuclei. 

Turning to ^^O, where an exact diagonalization of the 
18-body problem is out of reach, we compare our IM- 
SRG(2) results for the spectrum in Fig.|4]to calculations 
based on diagrammatic expansions (called the Q-box ex- 
pansion) commonly used to derive effective shell-model 
Hamiltonians [7[. In this context, one can also under- 
stand our choices for H^\ as follows. When one derives 
effective interactions among valence nucleons using per- 
turbation theory, then the many-body diagrams contain 
at least one vertex of F, 



pp'hh' •) J- pp' {yh or hv) •> 



orF 



{pq or qp)vv' • 



These interaction vertices are precisely the off-diagonal 
part driven to zero under the IM-SRG evolution. There- 
fore, the effective interactions among valence nucleons 
are directly given by PH f2{^P (o^ily at finite s, there 
would be perturbative corrections). For a clear compari- 
son to shell- model calculations for ^^O, we use the same 
empirical single-particle energies for the one-body part 
of Hq^^ in both IM-SRG and Q-box calculations, that is 
we replace the calculated one-body part in the IM-SRG 
by the empirical USDb [16] single-particle energies. 

The left panel of Fig. |4] compares the low- lying ^^O ex- 
citation energies obtained by diagonalizing the S(i-shell 
JT^eff derived from the IM-SRG (2) and the Q-box expan- 
sion. We also give the experimental energies [15], al- 
though good agreement with experiment is not required 
since three-nucleon (3N) forces are not included in the 
initial Hamiltonian and we do not fine-tune the single- 
particle basis to reproduce the experimental root-mean- 
square radius (see below) . All calculations are performed 
in the HF basis. For the perturbative Q-box results, the 
open symbols correspond to an i^eff that is calculated 
at first, second, and third-order, while the solid symbols 
include higher-order folded-diagram contributions to re- 
move the energy dependence of induced interaction ver- 
tices [7]. For the IM-SRG(2) results, as for ^Li, we ob- 
serve negligible differences in the calculated spectra for 
the two generators H^^ and i^2^. This implies that the 
truncation of the flow equations to two-body operators 
is a very good approximation. The IM-SRG (2) energies 
are similar to the results based on the perturbative Q- 
box expansion, where some differences from the "best" 
results (Q*^^^ plus folded- diagram contributions) are ex- 
pected because the IM-SRG is a nonperturbative method 
that includes many higher-order terms. 

The poor agreement with experiment in Fig. |4] com- 
pared to conventional shell-model calculations might be 
surprising. The reason for this discrepancy is that con- 
ventional calculations include additional phenomenology 
that improves agreement with experiment, but weakens 
the connection with the underlying Hamiltonian and mi- 
croscopic many-body theory. This can be understood 
from the right panel of Fig. HI which displays the hw de- 
pendence of the ^^O excitation energies for second-order 
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FIG. 4: Left panel: Excitation energies of ^^O obtained by diagonalizing the IM-SRG(2) ifeff in the S(i-shell, compared with 
results obtained at first (Q*^^-*), second (Cf^^^ and third (Q*^^-*) order in the Q-box expansion 0] for effective valence-shell 
interactions (the filled symbols include higher-order folded- diagram contributions). All results are for TujJ — 24 MeV, but a HF 
basis is used. For comparison, we also show the experimental energies [13. Right panel: Excitation energies of ^^O versus hjj 
calculated at the second-order Q-box level (plus folding) in a harmonic-oscillator and HF basis. For the results of both panels, 
an emax = 8 space was used, and to simplify the comparison, the calculations used single-particle energies from the USDb 
interaction [16'] for the diagonalization in the s(i-shell. 



Q-box calculations performed in the harmonic-oscillator 
(HO) and the HF basis. The HO-based spectrum ex- 
hibits a very strong huj dependence, while calculations 
in the HF basis are nearly independent of huj. Conven- 
tional calculations of i!/^eff work in the HO basis with 
hbj ^ 45 A~^/^ — 25 A~^/^ chosen to give the same root- 
mean-square radius as a sphere of uniform density. The 
fine-tuning of hu can therefore be understood as a phe- 
nomenological means to build in the correct saturation 
properties of nuclei, which are known to be deficient in 
ab-initio calculations starting from Hamiltonians with- 
out three- nucleon forces [9]. Our IM-SRG calculations 
present a microscopically-derived i^eff, but the incorrect 
saturation properties of the initial NN-only Hamiltonian 
translate into a poor description of the ^^O spectrum 
compared to empirical calculations carried out in a HO 
basis (at huj ^ 14 MeV for ^^O). 

Conclusions.- We have shown that the IM-SRG can 
be successfully generalized to open-shell systems and 
to a nonperturbative derivation of effective valence-shell 
Hamiltonians. The IM-SRG evolution decouples the 
physics of valence nucleons from the full Hilbert space, 
enabling exact diagonalizations in the valence space that 
are impossible in the full problem where all nucleons are 
active. First results were presented for ^Li, with ground- 
state and excited-state energies in very good agreement 
with ab-initio methods. We then applied the IM-SRG to 
^^O and compared our results to those obtained from con- 
ventional perturbative calculations of i^eff, demonstrat- 
ing that the IM-SRG provides a first viable nonperturba- 
tive approach to derive effective interactions for the shell 
model from nuclear forces. Work is in progress to ex- 
tend the IM-SRG to extended valence spaces, to effective 



operators and to the inclusion of three-nucleon forces. 
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